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Unitary ensembles of large n x n random matri- 
ces with a non-Gaussian probability distribution P [H] oc 
exp {—try [H]} are studied using a theory of polynomials or- 
thogonal with respect to exponential weights. Asymptotically 
exact expressions for density of levels, one- and two-point 
Green's functions are calculated. We show that in the large-n 
limit the properly rescaled local eigenvalue correlations are 
independent of P [H] while global smoothed connected corre- 
lations depend on P [H] only through the endpoints of spec- 
trum. We also establish previously unknown intimate con- 
nection between structure of Szego function entering strong 
polynomial asymptotics and mean-field equation by Dyson. 
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Random-matrix theory turned out to be general phe- 
nomenological approach to description of various phe- 
nomena in complex physical systems [0 such as com- 
plex nuclei, mesoscopic conductors or quantum chaotic 
systems. Due to complicated nature of these objects 
their treatment, based on the direct integration of mi- 
croscopic equations, runs across considerable obstacles. 
In this situation the more reasonable way is to deal with 
the joint distribution function P [H] of the matrix ele- 
ments of Hamiltonian H describing the system involved. 
If there is not preferential basis in the space of matrix 
elements (i.e. the system in question is "as ran- 
dom as possible" ) one has to require P [H] d [H] to be 
invariant under 7?,-transformation of the Hamiltonian: 
H 7?.^^H7?.. Here TZ is orthogonal, unitary or sym- 
plcctic matrix reflecting the fundamental symmetry of 
the underlying Hamiltonian. The general form of P [H] 
compatible with invariance requirement can be written 
as P [H] = exp {-trV [H]} with V [H] providing ex- 
istence of partition function Z. Transformation of H 
to diagonal form H = 5^^X5/3 and further integra- 
tion over orthogonal (/3 = 1), unitary (/3 = 2) or symplec- 
tic {P = 4) group dfj. (Sp) in the construction P [H] d [H] 
yield famous expression for the joint probability density 
function of the eigenvalues {x} of the matrix H: 



PiM) = ^ exp {-13 



i<3 



(1) 



Confinement potential V (x) contains, in an integral way. 



an information about correlations between different ma- 
trix elements of a random Hamiltonian H while logarith- 
mic interaction induced by the invariance requirement is 
usually treated as a pairwise repulsion between energy 
levels. We note that in Eq. (|l|) the symmetry parame- 
ter (3 is factored out from V {x) to fix density of levels 
in random-matrix ensembles with the same confinement 
potential but with different underlying symmetry. 

The structure of Eq. (Q) allows to represent exactly all 
the global and local statistical characteristics of the phys- 
ical system with broken time-reversal symmetry (/3 = 2) 
in the terms of polynomials orthogonal with respect to 
the measure da (x) = exp {—2V (x)) dx. 

Hitherto only few cases of strong confinement poten- 
tials were rigorously treated in the random-matrix the- 
ory. Random-matrix ensembles associated with classi- 
cal orthogonal polynomials have been considered in 
those related to some non-classical orthogonal polynomi- 
als have been treated in [Q. We also mention the 
basic work Q where ansatz for functional form of poly- 
nomial asymptotics was proposed to describe a class of 
non-Gaussian unitary random-matrix ensembles. 

Since generally an analytical calculation of the corre- 
sponding set of orthogonal polynomials is a non-trivial 
problem, a number of alternative approximate meth- 
ods were developed. The mean-field approximation pro- 
posed by Dyson allows to calculate density of levels in 
random-matrix ensemble. This approach being combined 
with the method of functional derivative of Becnakker 
makes it possible to compute global (smoothed) eigen- 
value correlations in large random matrices. Smoothed 
correlations can also be obtained by diagrammatic ap- 
proach of Brezin and Zee |^ or by invoking the argu- 
ments based on macroscopic electrostatics and linear re- 
sponse theory |^. We stress that all the methods men- 
tioned above allow to study correlations only in long- 
range regime and, in this sense, they are less informa- 
tive as compared with the method of orthogonal poly- 
nomials 1^]. It is worth pointing out the supersymmetry 
formalism ]To| , recently developed for invariant random- 
matrix ensembles with non-Gaussian probability distri- 
bution function P [H] , which is exceptional in that it al- 
lows to investigate local eigenvalue correlations and rep- 
resents a powerful alternative approach to the classical 
method of orthogonal polynomials. 

It is the aim of this Letter to outline how the prob- 
lem of non-Gaussian ensembles with unitary symmetry 
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can be handled by the method of orthogonal polynomi- 
als in rather general case. Extended version of the arti- 
cle containing detailed mathematical calculations will be 
published elsewhere . We stress that our treatment is 
exact (i.e. it does not involve any conjectures) and based 
on the recent results obtained in the theory of poly- 
nomials orthogonal with respect to exponential weights 
on R. 

Let us consider a class of even confinement potentials 
V (x) supported on the whole real axis x G (—00, +00) 
which are of smooth polynomial growth at infinity and in- 
crease there at least as |a;| ^^"^ {S is arbitrary small positive 
number). More precisely, we demand V (x) and (PV/dx^ 
be continuous in x G {0,+oo), and dV/dx > in the 
same domain of variable x. Further we define the func- 
tion T (x) = 1 + a; • d'^V/dx^ • {dV/dx)~^ and assume that 
for some A > 1 and B > 1 the inequality A < T (x) < B 
holds for X G (0, +00), and also for x positive and large 
enough x^ ■ \d^V/dx^\ ■ {dV/dx)" < const. The class 
of potentials V {x) satisfying all the above requirements 
is said to be of the Freud type. The typical examples of 
the Freud potentials are (i) V {x) = with a > 1; (ii) 
7 -I- X ) with a > 1, /3 e R, and 7 large 

enough. 

Defining the Freud (non-Gaussian) measure daj^ (x) = 
wjr (x) dx = exp {—2V (x)) dx we can introduce a set of 
orthogonal polynomials P„ (x) 



Pn (x) Pm {x) dajr (x) = S„ 



(2) 



for which the following basic result was obtained by D. 
S. Lubinsky 



lim 



-1-1 



dX 



[aU^Pn (a„A) 



1/2 



Re 



wjr (a„A) = 0. (3) 



Here parametrizations z = e and A = cos 9 were used. 

Szego function D {g; z), appeared in Eq. (||), is of fun- 
damental importance in the whole theory of orthogonal 
polynomials , and takes the form 



D {g; z) = cxp ( ^ 



1 -)_ ze"*''' 



1 - ze-^'fi 

The first argument of Szego function in Eq. (|^) is 

1/2 

Fn {(p) ^ exp {-V {an cos ip)) \sm ip\ ' , 



(4) 



(5) 



and a„ is the n-th Mhaskar-Rahmanov-Saff number being 
the positive root of the integral equation [H 



2a„ 



XdX 



V1-A2 \dxj ^^^^ 



dV 



(6) 



(In what follows it will be seen that a„ is none other 
than band edge for eigenvalues of corresponding random- 
matrix ensemble). 

Passing on to the integration over x = a„A in Eqs. 
(^) and (^, and leaving aside the problem of remainder 
term we arrive at the asymptotic formula, n ^ I, for 
orthogonal polynomials of the Freud type: 



Pnix) 



-Re 



z''D- 



F„ 



, X e (-a„, +an) . 

(7) 



Making use of the properties of Szego function D (g; e*^) 
|l5| one can obtain the following formula for the or- 
thonormal "wave functions" ipn {x) = Pn {x) exp{—V (x)) 
within the interval x € {--an, -f a„) in the large-n limit: 



Ipn (x) 



2 



1 - 



-1/4 



cos$„(x), (8) 



where 



*&n (x) — 7„ (x) 4- n arccos — 

\an 



In {X) - —V 
ZTT 



+a„ 



d^ 



vol 



h{0 



■In 



2F(0. 



(9) 



(10) 



(11) 



In Eq. ( p^ V stands for principal value of the inte- 
gral. Since the same formula takes place for the Erdos- 
type confinement potentials which grow faster than 
any polynomial at infinity (see |l7|] , Ch. 2) the further 
treatment applies equally to unitary random-matrix en- 
sembles with confinement potentials of the Freud- and 
Erdos-type. The simple examples of Erdos-type confine- 
ment potentials are (i) V (x) = exp^. (|x|") with a > 
and k > 1 (here exp^ denotes the exponent iterated k- 
times); (ii) V (x) — exp (in" (7 x^)) with a > 1, and 7 
large enough. 

All the global and local characteristics of the random- 
matrix ensembles are determined by two-point kernel |^ 



Kn (x, y) = 



kn-l i'n (y) i'n-l [x] - tpn (x) 1pn-l [v] 



y 



(12) 



where fc„ is a leading coefficient of the orthogonal poly- 
nomial Pn (x). Substituting Eq. (^ into Eq. ( p^ we 
obtain after straightforward calculations in the large-n 
limit: 
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Kn {x, y) = 



1 



TT{y -x) 



[fn{x)fn{y)]-'^' 



X \ COS $„ (x) COS $„ (y) 



X -y 



- sin $„ (y) cos $„ (x) V /" iv) 

+ sin $„ (x) cos $„ (y) v//„ (a;) } , (13) 

where (x) = 1 — (x/a„)^. This equation is vahd for 
arbitrary x and y lying within the band (— a„, a„). 

Equation ( |l3|) allows to compute smoothed (over the 
rapid oscillations) connected correlations i^c {x, y) P ] 

Vc{x,y) = {Vn {x)Un {y))-{Vn [x)) {Vn {v)) = -K^{x,y), 

where x ^ y and (...) denotes averaging over intervals 
I Ax I ttn and \Ay\ <^ a„ but still containing many 
eigenlevels. Direct calculations yield the simple univer- 
sal relationship 



Vc [x, y) 



xy 



[x-y) ^Jal - x'^^al - y' 



x^y 
(14) 



with dependence on the potential V (x) only through the 
endpoint a„ of the spectrum. 

Now we turn to the local properties of two-point kernel. 
Assuming that in Eq. ( |l3|) |a; — ?/| <C a„, and both x and 
y stay away from the band edge a„ we obtain 



Kn {x, y) 



sin [x) - jy)) 
7T{y-x) 



(15) 



This two-point kernel may be rewritten in locally univer- 
sal form noting that following representation for $„ (a;) 
takes place 

{x) = i arccos ( — ) - n [ uja„ (0 + T (2"- ^ 1) > 
2 \anj Jo 4 

(16) 

where 



l^a„ (x) 



2 ^ dV ^/^~~x^ ^^^^ 



e-^''di^f^—e 



Since the characteristic scale of the changing of (C) 
is ^ ~ {oj~ldLOa^ld^ ^ a„, Eq. (|l|) for |a; - y| < ? is 
reduced to universal form 



Kn {x, y) 



sin {i^Vn {v-x)\ 
Tr{y~x) 



(18) 



with 77„ — cj(j^ (^^T^) pl^-yi'^g the role of local density of 
levels. Correspondingly, local two-level cluster function 
being rewritten in rescaled variables s and s' 



Y2 {s, s') 



Kl{x,y) 



{Vn {x)) {Vn [y)) 



y=y( = ') 



siv? [tt (s — s' 



(19) 



proves universal Wigner-Dyson level statistics in the uni- 
tary random-matrix ensembles with Freud- and Erdos- 
type confinement potentials (here s ~ 77„x and s' = Vny 
are the eigenvalues measured in the local mean level- 
spacing) . 

Expression for density of levels {vn (x)) = Kn {x, x) im- 
mediately follows from Eq. ([T^), {vn (x)) = —-^d^n/dx. 
It can be written down in two equivalent forms. The first 



one, 



{vn {x = a„ cos 9)) 



TTGn Sm C 



d9 



n9 



SLTgD 



e |sin(p|;e'^ 



(20) 



establishes connection between density of levels in 
random-matrix ensembles with confinement potentials of 
Freud- and Erdos-type, and Szego function, Eq. (^, for 
corresponding set of orthogonal polynomials. 

Another representation of the density of states reads 



^d^ dV ^al - x"^ 



(21) 



Considering this expression as an equation for V (x) one 
can resolve it invoking the theory of integral equations 
with Cauchy kernel mM: 



Vix) 



dx' {vji [x')) In jx — x'l -I- 



(22) 



with /i being "chemical potential" . It is non more than 
famous mean-field equation which, in our treatment, fi- 
nally follows from the asymptotic formula Eq. (|^) for or- 
thogonal polynomials without any physical speculations. 
Thus, quite surprisingly, Szego function, Eq. (^, turns 
out to be closely related to the mean-field approximation 
by Dyson §]. 

Now the calculation of the one-point Green's function 



GP {x) 



tr- 



1 



= V 



(23) 



with p = ±1 becomes quite trivial. We obtain by means 
of Eqs. d^) and (H): 



dx TT .70 



idi dVy/^ 



e - x^ di 



(24) 
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Two-point connected Green's function 

where Xp = x + ipO and x^, ~ x' + ip'O {p,p' = ±1), can 
also be computed within the proposed formahsm. These 
rather cumbersome calculations lead to the following 
formula for the smoothed connected two-point Green's 
function: 

G^/ ix,x') = I \pp' , 



connected correlators and local universality of those with- 
out smoothing over rapid oscillations. The correlators 
were shown to depend on the measure only through the 
endpoints of spectrum (global universality) or through 
the local density of levels (local universality). Presented 
rigorous polynomial analysis enabled us to establish new 
local universal relationship in the random-matrix theory 
for normalized and rescaled connected two-point Green's 
function g^P (s, s') . It is also worthy of notice an interest- 
ing and quite surprising connection between the structure 
of Szego function and mean-field equation that has been 
revealed in the proposed formalism. In our opinion this 
fact is far from obvious. 

One of the authors (E. K.) gratefully acknowledges fi- 
nancial support of The Ministry of Science and The Arts 
of Israel. 



^\,x^x'. (26) 

It also can be shown that the two-point Green's 
function in the random-matrix theory exhibits the lo- 
cal universality. Introducing normalized and rescaled 
two-point connected Green's function g'^^ (s, s') = 
GPP {x,x') I {vn [x)) {vn [x' )) cxprcsscd in the terms of 
the eigenvalues s and s' measured in the local mean level- 
spacing, we obtain the following universal relationship 
that does not depend on the distribution function P [H] : 

+ i{p-p') !!!LllliZi^e*-(--')-g"(p-f ') . (27) 
(s - s') 

Note that expression of this type was previously obtained 
in [p^ only for Gaussian random-matrix ensemble using 
supersymmetry formalism. 

In summary, we have presented rigorous analytical con- 
sideration of the basic problems of the theory of large 
Hermitian random matrices with non-Gaussian distribu- 
tion function P [H] cx exp {— trF [H]}. Our treatment 
is applied to very wide class of monotonous confinement 
potentials which behave at least as \x\^^^ {5 > 0) and 
can grow as or even faster than any polynomial at in- 
finity (note, that this class of confinement potentials 
is much richer than that considered in ||]). We have 
shown that in such unitary random-matrix ensembles 
with non-Gaussian measure P [H] d [H] the density of lev- 
els and one-point Green's function essentially depend on 
the measure. In contrast, (connected) two-point char- 
acteristics of spectrum are rather universal. Indeed, we 
have observed global universality of smoothed two-point 
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